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Summary 

The  nature  of  Ag(lll)/Ti02  rutile  and  anatase  interfaees,  of  interest  for  the  design  of  memristors, 
has  been  studied  by  means  of  density  funetional  theory  (DFT)  ealeulations  using  various 
eomputational  approaehes.  We  have  eonsidered  interfaees  where  the  lattiee  mismateh  of  the  oxide 
erystalline  phase  and  the  metal  eleetrode  does  not  result  in  exeessive  strain.  The  bonding  at  the 
interfaee  is  very  weak  and  the  eharge  transfer  is  negligible  for  stoiehiometrie  oxides.  The  formation 
of  O  vaeaneies  has  a  lower  eost  at  the  interfaee  with  Ag  than  on  the  bare  titiania  surfaee  and  results 
in  stronger  adhesion  between  the  Ag  eleetrode  and  the  redueed  Ti02-x  oxide.  The  diffusion  of  Ag 
and  O  atoms/ions  aeross  the  interfaee  is  a  thermodynamieally  unfavorable  proeess  whieh  ean  oeeur 
only  at  high  temperatures  or  under  the  effeet  of  an  external  eleetrie  field.  Onee  Ag  atoms  are 
ineorporated  in  the  bulk  of  Ti02  they  ean  be  stabilized  in  interstitial  (more  favorable)  or 
substitutional  to  Ti  positions.  In  both  eases  Ag  is  ionized  and  transfers  the  valenee  eleetron  to  the 
host  erystal  with  formation  of  Ti^^  states.  The  Ag  atoms  remain  positively  eharged  even  when 
extended  Ag  ehains  are  formed  (nano filaments).  For  an  Ag  filament  inside  Ti02  to  exhibit 
eonduetive  behavior  a  higher  density  of  Ag  atoms  is  required  but  this  is  hardly  possible  in  the 
regular  bulk  erystalline  lattiee  of  Ti02  without  ereating  a  struetural  breakdown. 
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1. Introduction 

Five  main  technologies  are  presently  under  scrutiny  for  new  generation  of  non  volatile  memories: 
(1)  phase  change  memories  (PCM)  based  on  chalcogenide  glass/  (2)  memories  based  on  magnetic 
tunnel  junctions  and  the  magnetoresistive  effect  (MRAM),  ’  (3)  charge  trapping  memories  (CTM) 
where  information  is  stored  in  form  of  charge  traps  in  the  nitride  layer  of  a  semiconductor-nitride- 
semiconductor  capacitor  structure/  (4)  ferroelectric  memories  (FeRAM)/  and  (5)  transition  metal 
oxide  based  resistive  memories  (or  memristors,  RRAM)/"^^  An  aspect  common  to  all  these  devices 
is  that  a  material  is  sandwiched  between  two  electrodes  which  measure  the  resistivity  of  the 
material  hence  the  memory  state  of  the  device.  Changes  in  resistivity  can  be  generated  by  magnetic 
fields,  heat,  or  voltage  pulses,  switching  the  device  from  a  conductive  to  a  non  conductive  state. 

The  origin  of  the  resistive  switching  in  memristors  is  still  matter  of  debate  and  nice  reviews 
about  this  problem  can  be  found  in  ref  6  and  18.  It  is  assumed  that  typical  resistive  switching  is 
initiated  by  a  voltage-induced  partial  dielectric  breakdown  in  which  Joule  heating  modifies  the 
material  along  a  discharge  filament.  The  filament  may  be  composed  of  the  metal  electrode  which  is 
transported  into  the  insulator  or  due  to  the  formation  of  sub-oxides.  During  the  reset  transition  the 
conductive  filament  is  thermally  destroyed  by  the  high  local  power  density.  According  to  another 
mechanism  charge  injection  and  charge  displacement  effects  are  at  the  origin  of  the  phenomenon. 
The  idea  is  that  electronic  charges  are  injected  by  tunneling  at  high  fields  into  defect  sites  present  in 
the  insulator.  This  modifies  the  electrostatic  barrier  and  the  resistance  of  the  oxide.  In  a  variant  of 
this  model  charge  trapping  at  the  interface  can  affect  the  height  of  the  Schottky  barrier  at  a 
metal/semiconductor  interface.  A  third  model  is  that  charge  injection  may  result  in  a  insulator-to- 
metal  transition  by  effect  of  system  doping. 

Another  mechanism,  of  more  direct  interest  for  this  modeling  study,  is  the  possibility  that 
ionic  transport  and  red-ox  processes  are  at  the  basis  of  the  resistive  switching.  The  idea  is  that  the 
oxidation  of  the  active  metal  electrode  in  contact  with  the  oxide  phase  results  in  the  release  of  metal 
ions  (e.g.  Ag^).  This  is  followed  by  the  drift  of  the  mobile  ions  within  the  oxide  layer  and 
subsequent  discard  at  the  inert  metal  counter-electrode.  This  mechanism  should  lead  to  the 
formation  of  metal  dendrimers  inside  the  oxide  which  represent  the  conductive  filament  when  the 
state  of  the  cell  is  ON.  By  reversing  the  polarity,  the  filament  is  then  dissolved  resetting  the  system 
to  the  OFF  state.  A  variant  of  this  picture  is  that  what  diffuses  inside  the  insulating  phase  are  not 
cations  but  anions  moving  in  direction  of  the  anode  (actually,  what  is  migrating  are  probably 
oxygen  vacancies  towards  the  catode).  This  should  result  in  a  local  change  in  conductivity. 

Clearly,  there  are  many  open  questions  and  it  is  still  largely  unclear  which  type  of 
mechanism  is  responsible  for  the  observed  sudden  and  reversible  change  in  resistivity.  There  is  also 
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increasing  evidence  that  the  nature  of  the  interface  between  the  metal  electrodes,  typically  Ag  or  Pt, 
and  the  thin  oxide  dielectric  (NiO,  Ti02,  perovskites)  plays  an  important  role  in  the  process. 

Among  the  various  systems  proposed  there  are  Ti02  films  of  30-60  nm  thickness  deposited 
on  Pt  or  Ag  electrodes.  Titania  has  various  polimorphs,  with  rutile  being  the  thermodynamically 
most  stable  phase  while  anatase,  more  interesting  for  applications,  exists  as  a  metastable  phase 
below  about  800  °C.  However,  titania  films  can  also  be  amorphous,  or  can  exist  as  a  mixture  of 
amorphous  and  crystalline  phases,  with  the  amorphous  component  usually  increasing  by  decreasing 
the  film  thickness.  What  is  present  in  both  crystalline  and  amorphous  phases  is  the  short  range 
order,  as  shown  by  Raman  spectra.  Amorphous  films  are  preferred  in  order  to  eliminate  electrical 
leakage  along  grain  boundaries.  Other  requirements  for  ideal  materials  are  a  large  optical  band  gap 
and  thermodynamic  stability  at  high  temperatures. 

Despite  the  technological  importance  of  the  topic,  very  little  work  has  been  done  on  the  side 
of  electronic  structure  calculations.  In  a  recent  paper  the  rectifying  properties  of  Pt/Ti02  interfaces 
have  been  discussed  based  on  density  functional  theory  (DFT)  calculations,^^  but  there  are  not  many 
other  examples  in  the  literature.  In  this  work  we  have  considered  various  models  of  Ag/Ti02 
interfaces  using  a  first  principles  DFT  supercell  approach.  Since  Ti02  can  be  easily  reduced 
resulting  in  n-type  conductivity,  we  have  considered  both  stoichiometric  and  reduced  interfaces.  As 
it  will  be  shown  below,  these  latters  exhibit  much  larger  adhesion  energies  and  are  probably  more 
reactive.  We  also  considered  the  possible  interdiffusion  of  Ag  atoms  into  the  oxide  phase  or  of  O 
atoms  into  the  metal  electrode  in  order  to  verify  the  thermodynamic  cost  of  injecting  atomic  species 
at  the  metal/oxide  boundary.  In  a  second  part  of  the  study  we  have  considered  the  nature  of  Ag 
dopants  in  interstitial  and  substitutional  positions  of  bulk  Ti02  and  the  effect  of  increasing  the  Ag 
concentration  with  formation  of  Ag  nanofilaments  inside  the  crystal. 

2.  Methods,  assumptions,  procedures 

The  description  of  metal/oxide  interfaces  with  DFT  is  not  trivial  because  of  the  necessity  to  treat  at 
the  same  footing  both  the  metallic  (Ag)  and  the  insulating  (Ti02)  phases.  Furthermore,  the  self¬ 
interaction  problem  in  DFT  leads  to  an  underestimation  of  the  band  gap  of  Ti02  and  to  an  incorrect 
description  of  electron  localization  of  defects  like  oxygen  vacancies.^^  A  way  to  remove  the  self¬ 
interaction  problem  is  to  use  hybrid  exchange-correlation  functionals  where  the  Hartree-Fock 
exchange  is  mixed  in  various  amounts  with  the  DFT  exchange.  This  works  well  for  bare  oxides  but 
hybrid  functionals  do  not  describe  metals  very  welP,  limiting  their  use  in  metal/oxide  interfaces. 
Another  solution  is  represented  by  the  DFT+U  approach  where  a  set  of  atomic-like  orbitals  (here  Ti 
3d)  are  treated  with  a  new  Hamiltonian  which  depends  on  Ueff  =  U  -  J,  where  U  is  a  parameter 
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describing  the  energy  increase  for  an  extra  electron  on  a  particular  site  and  J  is  the  screened 
exchange  energy. 

For  this  work  we  used  a  multiple  strategy.  The  calculations  dealing  with  energetic  aspects  of 
the  Ag/TiOi  interfaces  have  been  done  at  the  generalized  gradient  approximation,  GGA  (PBE 
exchange-correlation  functional  ).  For  this  part  of  the  study  we  adopted  the  VASP  program  ’ 
which  uses  a  plane  wave  basis  set  and  a  projector  augmented  wave  method  (PAW)^^  for  the 
treatment  of  core  electrons.  A  TiOi  film  of  four  or  five  layers  thickness  was  interfaced  with  three 
Ag  layers;  further  details  are  reported  in  the  following  sections.  With  this  representation  of  the 
Ag(lll)  surface  both  the  lattice  constant,  4.17  A,  and  the  metal  work  function,  4.32  eV,  are 
reproduced  with  acceptable  errors  (4.09  A  and  4.46  eV  are  the  experimental  quantities, 
respectively).  The  atoms  within  the  supercell  have  been  relaxed  until  a  0.01  eV/A  convergence 
threshold  was  reached  for  atomic  forces.  Brillouin-zone  sampling  has  been  performed  on  a  4x1x1 
(rutile)  and  2x2x1  (anatase)  k-points  mesh  according  to  Monkhorst-Pack  integration  scheme.^^ 
Electronic  structure  properties  and  density  of  states  (DOS)  of  the  Ag/Ti02  interface  have  been 
obtained  in  selected  cases  using  a  DFT+U  approach  in  order  to  better  represent  the  band  gap  of 
Ti02  and  the  band  offset  of  the  metal  and  the  insulating  phases.  We  used  U  -  J  =  4  eV  since  this 
was  shown  for  Ti02  to  provide  results  similar  to  the  elassieal  B3LYP  hybrid  funetional  ’  used 
for  the  study  of  Ag  atoms  inside  the  bulk  of  the  oxide  (see  below).  With  this  U  value  the  band  gap 
of  bulk  anatase  Ti02  is  2.81  eV,  to  be  eompared  with  a  gap  of  2.25  eV  obtained  at  the  GGA  level 
and  of  3.2  eV  in  the  experiment. 

The  B3LYP  hybrid  funetional,  and  a  eompletely  different  eomputational  setup  were  used  to 
study  Ag  atoms  ineorporated  into  the  bulk  of  anatase  Ti02.  As  it  will  be  diseussed  below,  Ag  atoms 
aet  as  redueing  agents  with  formation  of  Ag^  ions  and  Ti^"^  states  whieh,  being  loealized,  give  rise  to 
eleetronie  levels  in  the  gap  of  Ti02,  an  aspeet  whieh  is  not  eorreetly  deseribed  by  pure  GGA 
approaehes.  We  performed  spin  polarized  B3LYP  ealeulations  using  the  CRYSTAL06  eode^^  sinee 
this  approaeh  has  been  sueeessfully  used  to  deseribe  redueed  Ti02.^^  The  Kohn-Sham  orbitals  have 
been  expanded  in  Gaussian  Type  Orbitals  (GTO)  and  the  all-eleetron  basis-sets  are  86411(d41)^^ 
for  Ti,  8411(dl)^'^  for  O,  and  311(d31)  for  Ag.^^  In  this  latter  ease  an  effeetive  eore  potential  was 
used.^^  We  eonsidered  two  different  models  of  bulk  anatase:  i)  a  nearly  eubie  2V2x2V2xl  supereell 
eontaining  32  Ti  and  64  O  atoms  and  ii)  a  3x1x2  supereell  eontaining  24  Ti  and  48  O  atoms.  The 
optimized  bulk  lattiee  unit  eell  parameters  were  taken  from  previous  B3LYP  ealeulations  {a  =  3.776 
A,  and  c  =  9.866  A).^^  The  band  gap  of  anatase  Ti02,  3.9  eV,  is  overestimated  eompared  to  the 
experimental  value,  3.2  eV.  Full  geometry  optimization  was  performed  until  the  largest  eomponent 
of  the  ionie  forees  was  less  than  SxlO'"^  a.u.  The  reeiproeal  spaee  is  sampled  aeeording  to  a  regular 
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sublattice  with  a  shrinking  factor  (input  IS)  of  2  for  the  larger  supercell  and  of  4  for  the  smaller  one. 
The  DOS  curves  have  been  computed  with  a  4  x  4  x  4  k-point  mesh. 

3.  Results  and  discussion 

In  real  devices  the  interface  between  the  metal  electrode  and  the  insulating  film  involves  a 
poly  crystalline  material  (Ag)  and  an  amorphous  phase  (Ti02).  In  our  modeling  study,  in  order  to 
reduee  the  complexity  of  the  system,  we  consider  metal  and  oxide  single  crystal  surfaces.  While  this 
allows  one  to  have  a  microscopic  view  of  the  interface  structure,  it  may  also  result  in  other 
problems  like  the  need  to  create  commensurate  interfaces  where  little  strain  remains  after  adapting 
the  Ag  lattice  to  that  of  titania. 

3.1  Ag/Ti02  intefface  structures 

There  are  not  too  many  studies  dealing  with  stoichiometric  epitaxial  films  of  rutile  and  anatase  Ti02 
on  metal  single  crystal  surfaces.  On  the  contrary,  intensive  work  has  been  reported  by  the  group  of 
Granozzi  on  new  phases  associated  to  ultrathin  titania  films  on  Pt.^*  Reeently,  a  monolayer  Ti02 
film  with  lepidocrocite  structure  grown  on  Ag(lOO)  has  been  reported.^^  On  the  theoretical  side, 
Lopez  and  Norskov"^®  studied  the  Au(l  1 1)/Ti02(l  10)  rutile  interfaee  which  however  implies  a  large 
modification  of  the  lattiee  parameter  of  the  gold  metal  in  order  to  create  a  commensurate  interface. 
This  configuration  has  been  used  as  a  starting  point  for  our  study  of  the  Ag(lll)/TiO2(110)  rutile 
interface;  the  Ag(l  1 1)  slab  has  been  compressed  by  7%  in  the  [1-1 1]  direction  and  expanded  by  the 
same  amount  in  the  [001]  direction.  The  supercell  used  is  2x1  and  contains  18  Ag  atoms  and  16 
Ti02  units.  The  Ag  atoms  have  been  placed  either  on-top  of  the  bridging  oxygens  of  rutile  (110), 
Fig.  la,  or  in  the  regions  between  them.  However,  due  to  the  large  lattice  mismatch,  when  the 
structure  is  relaxed  unphysical  distortions  occur  in  the  metal  phase.  In  particular,  a  substantial 
rumpling  of  the  Ag(lll)  planes  occurs  due  to  the  strain  in  the  slab.  This  is  shown  by  the  fact  that 
the  same  distortion  is  found  also  for  an  unsupported  Ag(lll)  film  with  the  same  lattice  constants. 
Due  to  the  substantial  strain  in  the  structure,  this  model  will  be  no  longer  discussed. 
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Fig.  1  -  (a)  The  Ag(l  1 1)/Ti02(l  10)  rutile  interface  (0-top  configuration),  (b)  The 
Ag(l  1 1)/Ti02(100)  rutile  interface  (0-top  configuration). 

To  reduce  the  strain  we  have  considered  the  next  most  stable  rutile  surface,  the  (100)  one.  Flere  one 
can  reduce  the  compression  of  the  Ag  film  along  the  [1-11]  direction  below  5%.  The  expansion 
along  the  [001]  direction  remains  of  7%.  The  supercell  used  is  3x1  and  contains  18  Ag  atoms  and 
15  Ti02  units.  Also  in  this  case  we  considered  a  starting  configuration  where  the  Ag  atoms  are 
above  the  bridging  oxygens  of  rutile  (100),  Fig.  lb. 

During  geometry  relaxation,  the  metal  film  undergoes  a  large  structural  rearrangement  with 
distortions  along  the  z  direction  of  the  Ag  planes.  However,  differently  from  the  previous  case,  this 
distortion  is  due  to  the  chemical  interaction  of  the  metal  with  the  oxygen  atoms  of  the  TiOi  surface; 
the  standing  along  Ag  film  does  not  exhibit  the  same  kind  of  reconstruction.  The  adhesion  energy, 
Eadh,  defined  as  the  energy  of  the  interface  structure  minus  that  of  the  separate  components  in  their 
optiman  geometry,  is  low,  25.0  meV/A^,  the  charge  transfer  from  Ag  to  Ti02  is  negligible,  9.8  •  10'^ 
e/A^,  and  the  interface  distance,  defined  as  difference  of  the  average  vertical  positions  of  the  Ag  and 
of  the  O  or  Ti  top  layers,  2.89  A,  is  rather  long,  Table  1.  As  shown  by  the  DOS  curves  obtained  at 
the  DFT-i-U  level.  Fig.  2a,  there  is  almost  no  hybridization  of  the  Ag  and  TiOi  states,  which  is 
consistent  with  the  very  small  charge  transfer  at  the  interface.  The  Fermi  level  lies  high  in  the  band 
gap  of  the  oxide.  The  interface  represents  a  Shottky  contact  with  a  conduction  band  offset  from  the 
Fermi  energy  of  0.8??  eV.  However,  due  to  the  extremely  weak  interaction  and  the  absence  of  O- 
Ag  hybridization,  there  are  almost  no  metal  induced  gap  states  (MIGS). 


Fig.  2  -  DOS  curves  of  (a)  the  stoichiometric  Ag(l  1  l)/TiO2(100)  rutile  interface;  (b)  the  reduced 
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Ag(l  1 1)/Ti02(100)  rutile  interface  (DFT+U  results) 


It  is  well  known  that  titania  is  easily  reduced.  For  this  reason  we  have  considered  also  the  nature  of 
the  interface  between  a  highly  reduced  rutile  (100)  surface  and  Ag(lll).  All  the  bridging  oxygens 
of  rutile  (100)  have  been  removed  from  the  surface  and  the  properties  of  the  Ag(lll)/TiO2-n(100) 
rutile  interface  have  been  recalculated.  The  energy  required  to  removed  oxygen  from  the  Ti02 
phase  and  form  a  reduced  Ag/Ti02-n  interface,  Fred,  normalized  to  the  number  of  removed  O  atoms 
per  unit  cell,  n,  is  defined  as: 

Ered(Ag/Ti02)  =  [E(Ag/Ti02-n)  +  n/2  E(02)  -  E(Ag/Ti02)]/n  (1) 

It  should  be  compared  with  the  same  process  occurring  on  the  bare  TiO2(100)  surface: 

Ered(Ti02)  =  [E(Ti02.„)  +  n/2  E(02)  -  E(Ti02)]/n  (2) 

While  Ered(Ti02)  is  about  4.3  eV  per  O  atom  (small  fluctuations  are  connected  to  the  size  of  the 
supercell  and  the  number  of  Ti02  layers),  Ered(Ag/Ti02)  is  3.84  eV  only,  indicating  that  the 
formation  of  a  metal/oxide  interface  favors  the  chemical  reduction  of  the  oxide.  Furthermore,  it  has 
been  shown  in  the  literature  and  confirmed  in  this  study,  that  the  oxygen  removal  is  easier  for  the 
titania  layer  in  direct  contact  with  the  metal  than  for  internal  bulk  layers.  This  result  is  quite  general 
and  has  been  demonstrated  for  other  metals  and  oxides.^®’"^* 


Table  1  -  Properties  of  Ag(l  1 1)/Ti02  interfaces  (Eadh  =  adhesion  energy,  CT  =  charge  transfer  from 
Ag  to  Ti02,  Zinterface  =  interface  distance^) 


supercell 

Eadh,  meV/A^ 

CT,  e/A^ 

^interface,  A 

Ag(lll)/Ti02 

Rutile  (100) 

3x1 

25.0 

bo 

0 

(A) 

2.89 

Ag(lll)/Ti02-x 

Rutile  (100) 

3x1 

62.3 

-29.4  •  10'^ 

2.49 

AgOx(lll)/Ti02-x 

Rutile  (100) 

3x1 

84.0 

- 

2.48 

Ag(lll)/Ti02 

Anatase  (001) 

2V2  X  2V2  r45° 

13.1 

9.5  •  10“^ 

2.59 

a)  defined  as  difference  of  the  average  vertical  positions  of  the  Ag  and  of  the  O  or  Ti  top  layers 


Not  surprisingly,  the  Ag/Ti02  interaction  becomes  much  stronger  for  the  reduced  surface. 
The  adhesion  energy,  62.3  meV/A^,  is  more  than  twice  that  for  the  stoichiometetric  interface,  the 
charge  transfer  is  more  pronounced  and  goes  in  the  opposite  direction,  from  Ti02  to  Ag,  -29.4  •  10'^ 
e/A^,  and  the  interface  distance  is  shorter,  2.49  A,  Table  1.  This  has  direct  consequences  on  the 
DOS  curves  and  in  particular  results  in  the  presence  of  MIGS  in  the  band  gap  of  the  oxide  which 
are  now  more  pronounced.  Fig.  2b.  The  Fermi  level  is  now  above  the  tails  of  the  bottom  of  the 
conduction  band,  indicating  a  partial  metallization  of  the  rutile  (100)  surface  at  this  level  of 
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chemical  reduction. 

So  far  we  have  assumed  that  during  the  reduction  of  Ti02  the  O  atoms  recombine  and  leave 
the  sample  in  form  of  O2  molecules,  see  eq.  (1)  and  (2).  However,  at  a  real  interface  it  is  possible 
that  the  O  atoms  are  displaced  from  their  lattice  positions  and  diffuse  into  the  bulk  of  the  metallic 
phase.  This  process  could  lead  to  an  intermediate  AgxOy  phase  in  contact  with  a  reduced  titania 
layer.  A  process  of  this  kind  could  facilitate  the  generation  of  Ag^  ions  to  be  injected  into  the  oxide 
phase.  To  this  end,  we  have  considered  a  simple  model  where  the  bridging  O  atoms  of  the 
TiO2(100)  rutile  surface  have  been  removed  and  incorporated  into  the  Ag  metal  in  interstitial 
positions  between  the  first  and  the  second  layers  of  Ag(l  1 1).  We  have  considered  two  supercells  of 
different  dimensions,  a  3x1  cell,  Fig  3a,  and  a  3x2  supercell,  Fig.  3b.  In  the  first  case  all  the 
bridging  O  atoms  were  moved  into  the  Ag  metal,  while  in  the  second  case  only  one  every  two  O 
atoms  has  been  displaced.  Fig.  3b. 


k  " 


Fig.  3  -  Formation  of  a  reduced  Ag(l  1 1)/Ti02(100)  rutile  interface  and  of  an  oxygen  layer  inside 

the  Ag  electrode 

This  process  leads  to  a  partly  reduced  titania  phase,  Ti02-n,  and  a  weakly  oxidized  metal  phase, 
AgOn.  The  corresponding  formation  energy,  Ered-ox,  normalized  to  the  number  of  displaced  O  atoms, 
is  defined  as: 


Ered-ox(Ag/Ti02)  =  [E(Ag0n/Ti02-„)  -  E(Ag/Ti02)]/n  (3) 

After  relaxation,  the  O  atoms,  initially  placed  in  octahedral  cavities,  have  moved  into  (distorted) 
tetrahedral  cavities.  This  is  also  the  site  occupied  by  O  atoms  in  the  native  Ag20  oxide.  Ered-ox  is 
still  positive  and  rather  large.  For  the  small  supercell,  corresponding  to  high  O  concentration,  Ered-ox 
is  3.73  eV;  for  comparison,  Ered  value  for  the  Ag/Ti02  interface,  eq.  (1),  is  3.84  eV,  i.e.  only  slightly 
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larger  indicating  little  gain  when  O  atoms  are  injected  into  the  Ag  electrode  compared  to  the 
formation  of  gas-phase  O2.  For  the  larger  supercell  (smaller  concentration  of  O  interstitials),  Ered-ox 
decreases  to  3.03  eV,  still  indicating  an  endothermic  process.  The  reason  can  be  understood  by 
considering  the  low  affinity  of  Ag  for  O.  If  we  place  an  O  atom  in  the  bulk  of  Ag  and  we  take  V2O2 
as  a  reference,  we  find  that  the  process  is  endothermic  by  0.47  eV  and  0.65  eV  for  an  O  atom 
occupying  an  octahedral  and  a  tetrahedral  cavity,  respectively.  Clearly,  it  is  not  possible  to 
compensate  the  high  energy  cost  to  remove  oxygen  from  the  titania  surface  by  incorporating  it  into 
silver.  The  consequence,  however,  is  that  the  adhesion  energy,  computed  as  the  energy  required  to 
separate  the  AgOn  and  the  Ti02-n  layers,  increases  to  86.4  meV/A^,  and  that  an  AgOn  phase  is 
formed.  Here  the  Ag  atoms  are  partly  oxidized,  an  effect  which  could  favor  the  generation  of  Ag^ 
ions  and  their  diffusion  into  the  oxide. 

In  order  to  study  these  processes  we  have  considered  another  titania  surface,  the  (001)  face 
of  anatase,  which  has  the  characteristic  to  be  nearly  commensurate  with  the  Ag(lll)  surface.  The 
(001)  is  the  second  most  stable  anatase  surface,  and  it  has  been  chosen  because  one  can  identify  a 
square  unit  cell  where  the  lattice  mismatch  at  the  interface  with  Ag(l  1 1)  is  of  1.8%  only.  This  2^2  x 
2V2  r45°  unit  cell  contains  141  atoms  (45  Ag  and  32  Ti02  units).  Fig.  4a-c,  and  is  about  twice  that 
used  for  the  Ag(l  1  l)/rutile  interfaces. 
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Fig.  4  -  The  Ag(l  1  l)/TiO2(001)  anatase  interface,  (a)  Surface  unit  cell  of  anatase  (001);  (b)  surface 
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unit  cell  of  Ag(l  1 1);  (c)  side  view;  (d)  Diffusion  of  an  O  atom  from  the  oxide  to  the  metal  phase. 


After  optimizing  the  geometry  there  is  no  sign  of  rumpling  or  distortions  inside  the  Ag  layer.  This  is 
also  because  the  interaction  of  the  two  components  is  quite  weak,  with  an  adhesion  energy  of  13.1 
meV/A^  only,  a  charge  transfer  from  Ag  to  Ti02  of  9.5  •  10'^  e/A^,  and  an  interface  distance  of  2.59 
A,  Table  1.  The  DOS  eurves  obtained  at  the  DFT+U  level  are  given  in  Fig.  5. 


ElaV) 


Fig.  5  -  DOS  Ag(l  1 1)/Ti02  anatase  (DFT+U  result) 

As  for  the  stoichiometric  Ag(l  1 1)/Ti02(100)  rutile  interface,  there  is  almost  no  hybridization  of  the 
Ag  and  Ti02  states,  which  is  consistent  with  the  negligible  charge  transfer.  The  Fermi  level  lies 
about  2  eV  above  the  top  of  the  O  2p  valence  band.  However,  due  to  the  extremely  weak 
interaction,  the  are  no  metal  induced  gap  states.  In  short,  this  interface  structure  has  virtually  no 
strain  but  has  also  no  chemical  interactions  at  the  interface,  making  injections  of  atoms  or  ions  into 
one  of  the  two  phases  more  difficult. 

This  aspect  has  been  investigated  by  studying  the  same  processes  described  above  for  the 
Ag-rutile  surface.  In  this  case,  however,  thanks  to  the  large  dimensions  of  the  supercell,  we  can 
model  the  displacement  of  isolated  O  atoms  of  the  Ti02  surface  and  their  incorporation  into  the  Ag 
electrode.  The  first  step  has  been  the  calculation  of  the  formation  energy  of  an  oxygen  vacancy  on 
the  bulk  terminated  TiO2(001)  anatase  surface.  This  is  2.93  eV  when  the  O  atom  is  removed  from 
the  first  layer  and  3.98  eV  when  O  is  taken  from  a  bulk  layer  (notice  that  the  (001)  anatase  surface 
in  vacuum  undergoes  a  strong  relaxation  so  that  the  formation  energy  of  a  vacancy  for  the  relaxed 
surface  is  4.3-4. 6  eV,"^^  depending  on  the  layer  that  is  reduced).  The  low  value  of  the  vacaney 
formation  energy  for  the  top  layer  of  the  unreconstructed  TiO2(001)  surface  reported  above,  2.93 
eV,  is  thus  somewhat  artificial  and  a  better  estimate  is  that  computed  when  oxygen  is  removed  from 
subsurface  layers,  nearly  4  eV. 

We  considered  the  Ag(lll)/TiO2(001)  interface  and  studied  the  cost  of  reducing  titania  and 
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simultaneously  oxidizing  silver.  This  was  done  by  displacing  a  single  O  atom  from  TiOi  and 
including  it  in  an  interstitial  subsurface  site  of  Ag(lll),  Fig.  4d.  The  total  energy  of  this  rex-ox 
system  is  2.72  eV  higher  than  the  non-defective  interface.  Depending  on  which  O  atom  is  removed, 
this  energy  changes  slightly  but  never  more  than  0.1  eV.  Thus,  also  in  this  case  we  find  that  the 
process  is  endothermic  and  can  occur  only  at  high  temperatures  or  under  the  effect  of  an  high 
voltage  so  that  substantial  energy  barriers  can  be  overcome.  On  the  other  hand,  the  energy  required 
to  move  O  into  the  metal  is  lower  than  that  needed  to  form  14  O2  in  gas-phase.  Therefore,  any 
reduction  process  occurring  once  the  interface  is  formed  will  result  in  partial  oxidation  of  the  Ag 
layer  rather  than  in  release  of  gas-phase  oxygen. 

When  an  O  vacancy  is  created  on  anatase  TiO2(001),  an  Ag  atom  is  strongly  displaced  from 
the  interface  towards  the  vacancy.  Fig.  4d.  This  Ag  atom  is  slightly  positively  charged,  according  to 
the  Bader  analysis,  suggesting  that  the  presence  of  O  vacancies  at  the  interface  could  represent  a 
channel  for  the  diffusion  of  positive  Ag  ions  into  the  oxide. 

In  a  similar  way,  we  considered  the  energy  required  to  move  an  Ag  atom  from  the  metal  into 
the  oxide,  i.e.  the  inverse  process  with  respect  to  the  O  diffusion.  To  this  end,  an  Ag  atom  taken 
from  the  perfect  interface  has  been  included  in  a  interstitial  cavity  of  anatase.  We  found  that  the  Ag 
species  is  ionized,  as  shown  by  a  Bader  charge  of  +0.52  (see  also  §  3.2).  However,  during  the 
geometry  optimization  the  Ag  atom  moves  back  to  its  original  position  restoring  the  non-defective 
interface.  This  is  a  clear  sign  that  the  process  is  rather  endothermic,  and  that  diffusion  of  Ag  is 
possible  only  under  the  effect  of  external  driving  forces  like  those  provided  by  an  applied  voltage. 

3.2  Nature  of  Ag  atoms  incorporated  into  bulk  anatase 

The  presence  of  Ag  atoms  inside  the  bulk  of  anatase  Ti02  and  their  effect  on  the  electronic  structure 
of  the  system  has  been  investigated  by  performing  supercell  calculations  with  the  CRYSTA06  code 
and  the  B3LYP  exchange-correlation  functional.  We  considered  two  possible  sites  where  the  Ag 
atoms  are  placed,  substitutional  to  Ti  or  in  interstitial  cavities.  When  Ag  is  substitutional  to  Ti  it  is 
coordinated  to  six  neighboring  oxygens  via  two  long  apical  and  four  short  equatorial  bonds  of  2.22 
and  2.03  A  respectively.  The  local  structure  remains  that  of  undoped  anatase  except  for  the  bond 
lengths  that  are  about  0.2  and  0. 1  A  longer.  While  Ti  in  Ti02  is  formally  +IV,  the  Ag  atom  can  only 
partly  supply  the  electrons  needed  to  fill  the  valence  of  the  O  atoms.  When  Ag  substitutes  a  Ti  atom, 
Fig.  6a,  is  formally  in  a  +II  oxidation  state  (3d^  configuration),  with  an  hole  in  the  3dz^  orbital.  The 
fact  that  Ag  is  only  a  two-electron  donor,  results  in  two  holes  in  the  2p  levels  of  the  O  atoms,  one 
fully  localized  and  one  more  delocalized.  Fig.  6b.  Both  doublet  and  quartet  spin  configurations  have 
been  considered,  with  this  latter  only  slightly  lower  in  energy  (by  0.03  eV).  However,  the  energy 
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cost  to  remove  a  Ti  atom  from  the  lattice  and  bring  it  to  the  gas-phase  and  replace  it  with  Ag  is 
extremely  high,  22.6  eV.  Of  course,  a  much  lower  energy  is  expected  if  the  Ti  atom,  instead  of  being 
removed  from  the  crystal,  is  displaced  into  an  interstitial  cavity. 


Fig.  6  -  (a)  Structure  and  (b)  spin  distribution  of  a  Ag  atom  substitutional  to  Ti  in  bulk  anatase;  (c) 
structure  and  (d)  spin  distribution  of  an  interstitial  Ag  atom  in  bulk  anatase 


There  are  various  possible  positions  where  the  Ti  atom  can  be  accommodated  and  their  total 
energies  have  been  tested.  In  the  most  stable  site,  Ti  is  coordinated  to  five  O  atoms.  The  overall 
energy  cost  to  replace  a  lattice  Ti  atom  with  an  Ag  atom  and  to  displace  it  into  an  interstitial  site  is 
2.4  eV,  an  energy  which  can  be  easily  provided  during  thermal  treatment  of  the  sample  or  under  an 
applied  voltage. 

Quite  interesting  from  the  point  of  view  of  the  stability  is  the  Ag  interstitial  species,  Fig.  6c. 
The  formation  energy,  defined  as  the  energy  required  to  include  a  gas-phase  Ag  atom  into  the  bulk 
of  anatase,  is  1.1  eV,  indicating  a  moderately  endothermic  process.  However,  the  structure  where  Ag 
is  interstitial  and  Ti  occupies  a  lattice  position  is  1.3  eV  lower  in  energy  than  that  described  above 
where  the  two  atoms  have  interchanged  their  positions.  This  shows  unambiguously  the  preference 
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of  Ag  for  interstitial  sites.  In  the  optimized  strueture  Ag  assumes  a  distorted  square  planar 
eoordination  with  O  atoms  at  about  2.2  A  while  the  apieal  O  atom  is  2.93  A  apart,  Fig.  6e.  From  the 
eleetronie  point  of  view,  there  is  elear  evidenee  that  the  Ag  atom  aets  as  a  redueing  agent  forming  a 
Ag^  ion.  Fig.  7.  The  Ag  5s  eleetron  is  transferred  to  the  bottom  of  the  eonduetion  band  whieh,  being 
formed  of  Ti  3d  states,  leads  to  the  formation  of  a  loealized  Ti^^  (3dxy)^  speeies.  Fig.  6d  and  Fig.  7. 
In  partieular,  the  eleetron  is  loealized  on  a  Ti  ion  nearest  neighbor  of  the  interstitial  Ag. 


ENERGY  (eV) 

Fig.  7  -  Density  of  states  of  Ag  interstitial 


These  results  provide  elear  evidenee  that  an  isolated  Ag  atom  inside  anatase  Ti02  prefers 
interstitial  positions,  looses  the  5s  valenee  eleetron,  and  is  in  a  positive  oxidation  state.  Only  if  the 
density  of  the  ineorporated  atoms  inereases,  they  may  switeh  to  a  neutral  state  whieh  eould  lead  to 
the  formation  of  a  metallie  “filament”.  This  meehanism  will  eompete  with  the  tendeney  of  Ti  3d 
orbitals  to  trap  eleetrons  and  form  Ti^^  states. 

In  order  to  verify  this  possibility  we  have  designed  a  different  3x1x2  supereell  (see 
Computational  method)  whieh  allows  one  to  have  Ag  atoms  at  shorter  distanees  (the  Ag-Ag 
distanees  are  3.78  A).  This  eould  result  in  the  formation  of  a  band  due  to  the  overlap  of  the  Ag  5s 
levels  and,  possibly,  in  the  formation  of  a  metallie  wire.  We  performed  this  test  for  both 
substitutional  and  interstitial  Ag,  Fig.  8.  For  the  filament  model  ereated  by  the  repetition  of 
substitutional  Ag  atoms  we  found  a  similar  eleetronie  and  geometrie  strueture  as  for  an  isolated 
speeies.  Ag  has  a  distorted  oetahedral  eoordination  and  the  formation  energy  is  still  very  high.  The 
absenee  of  a  direet  orbital  interaetion  among  the  Ag  atoms  indieates  that  it  is  not  possible  to  form  an 
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Ag  nanowire  by  substitutional  defects  mainly  because  the  Ag-Ag  separation  is  still  too  large. 


Fig.  8  -  (a)  Chain  of  substitutional  Ag  atoms  in  anatase  Ti02  and  (b)  chain  of  interstitial  Ag  atoms 

in  anatase  Ti02  (the  plot  shows  the  spin  density). 


Also  the  chain  of  interstitial  Ag  atoms  does  not  exhibit  a  metallic  character.  The  nature  of  the 
individual  components  is  that  of  an  Ag^  ion.  The  only  differences  with  respect  to  the  isolated  Ag 
case  are  that  (i)  from  a  structural  point  of  view  there  is  a  displacement  at  the  center  of  square  planar 
pyramid  of  the  Ag  atom  that  becomes  five-coordinated  and  (ii)  from  the  side  of  the  electronic 
structure  the  unpaired  electron  is  localized  on  two  Ti  atoms  and  occupies  a  3dxy  atomic  orbital  of 
one  Ti  neighbor  and  a  3dyz  orbital  of  the  other  Ti  neighbor,  Fig.  8b.  With  respect  to  an  isolated  Ag 
atom  the  formation  of  the  chain  leads  to  a  small  energy  gain  of  0.39  eV/Ag  atom,  suggesting  that  a 
small  but  non  negligible  attractive  interaction  occurs  among  the  Ag  interstitials  and  that  the 
inclusion  of  more  silver  could  result  in  the  formation  of  a  metallic  wire. 

The  amount  of  Ag  incorporated  in  the  oxide  has  been  increased  by  including  two  Ag  atoms 
in  the  same  supercell,  one  substitutional  to  Ti,  Agsub,  and  one  in  interstitial  positions,  Agint,  Fig.  9.  In 
this  way  the  Ag-Ag  distance  becomes  2.75  A,  not  far  from  that  of  bulk  silver,  2.88  A.  The  resulting 
Ag  wire  has  a  zig-zag  structure,  Fig.  9.  It  is  interesting  to  note  that  the  same  structure,  i.e.  a  zig-zag 
chain,  is  adopted  by  a  free  Ag  wire,  not  embedded  in  the  oxide.  This  structure  is  more  favorable 
than  a  linear  Ag  chain;  furthermore,  the  geometry  optimization  leads  to  small  changes  in  the  Ag-Ag 
distances  (2.61  and  2.65  A)  which  are  typical  of  a  Peierls  distortion  in  one-dimensional  conductors. 
The  free  Ag  zig-zag  wire  exhibits  5s  density  of  states  at  the  Fermi  level,  typical  of  a  metal.  Things 
are  quite  different  when  we  consider  the  Ag  wire  inside  Ti02.  In  fact,  in  the  DOS  (not  reported)  we 
observe  the  appearance  of  an  occupied  state  (spin  polarized,  a  component)  with  O  character  just 
above  the  top  of  the  Ti02  valence  band  and  close  to  the  Fermi  level.  This  is  because  the  Ag  atoms 
are  in  positive  oxidation  state,  each  acting  as  one-electron  donor.  However,  the  two  Ag  atoms  in  the 
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supercell  are  not  sufficient  to  compensate  for  the  missing  Ti  atom  (one  of  the  two  Ag  atoms  replaces 
a  lattice  Ti  atom).  As  a  consequence,  two  holes  are  present  in  the  O  2p  states,  partly  mixed  with  the 
Ag  4d  states.  Thus,  despite  the  formation  of  a  continuous  Ag  wire  there  is  no  sign  of  the  appearance 
of  associated  conductive  states.  The  Ag  5s  band  remains  empty  and  well  above  the  bottom  of  the 
Ti02  conduction  band. 


Fig.  9  -  Zig-zag  chain  of  Ag  interstitial  and  Ag  substitutional  in  bulk  anatase.  (a)  front  view;  (b) 

side  view. 

To  remove  the  problem  of  electron  deficiency  originating  from  the  replacement  of  one  Ti 
with  two  Ag  atoms  per  unit  cell,  we  have  considered  a  structure  where  also  an  O  vacancy  has  been 
created.  The  rational  is  that  the  electron  deficiency  (holes  in  the  valence  band  of  Ti02)  can  be 
compensated  by  creating  an  O  vacancy,  a  two-electron  center  associated  to  the  formation  of  two 
Ti  ions.  Indeed,  when  we  create  the  O  vacancy  there  are  no  more  holes  in  the  Ti02  valence  band 
and  the  system  is  electrically  compensated.  This  is  relevant  because  we  have  now  a  system  where  a 
continuous  chain  of  Ag  atoms  has  been  created  without  perturbing  in  a  major  way  the  electronic 
structure  of  the  anatase  host  crystal.  If  one  now  injects  one  electron  in  the  system,  there  are  two 
possibilities:  (1)  the  extra  electron  enters  in  the  band  generated  by  the  Ag  5s  levels,  thus  giving 
metallic  character  to  the  Ag  wire;  (2)  the  electron  prefers  to  occupy  the  Ti  empty  3d  states  which 
contribute  to  the  bottom  of  the  conduction  band  of  Ti02.  This  possibility  has  been  checked  by 
performing  a  calculation  for  the  Ag-doped  anatase  supercell  where  one  extra  electron  has  been 
added  (this  results  in  a  charged  supercell  which  requires  the  use  of  a  background  of  charge  in  the 
calculation).  The  analysis  of  the  DOS  clearly  shows  that  the  extra  electron  prefers  to  occupy  Ti  3d 
states  (formation  of  Ti^^)  instead  of  the  Ag  5  s  band  which  remains  higher  in  energy. 
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In  conclusion,  there  is  a  strong  tendency  for  silver  to  be  stabilized  as  positive  ion  inside  the 
bulk  of  anatase,  and  this  remains  true  even  for  densely  packed  one-dimensional  chains.  This  means 
that  the  screening  of  the  positive  charge  from  the  surrounding  oxide  is  rather  efficient  so  as  to  allow 
the  presence  of  an  array  of  ions  in  this  configuration.  Probably  a  higher  concentration  of  Ag  inside 
the  material  is  needed  in  order  to  reach  metallization. 

4.  Conclusions 

In  this  work  we  have  considered  simplified  models  of  Ag(lll)/Ti02  interfaces.  These  systems  are 
of  relevance  for  the  design  of  memristor  devices.  The  focus  has  been  on  the  following  aspects. 

(1)  We  first  considered  models  of  the  Ag/Ti02  interface  where  the  lattice  mismatch  between 
the  two  phases  is  below  a  few  %  in  order  to  reduce  the  strain.  This  is  important  in  our  calculations 
since  we  are  dealing  with  single  crystal  phases  and  the  presence  of  a  large  strain  in  the  material  can 
result  in  unphysical  structural  distortions.  We  have  found  that  the  anatase  (001)  face  has  a  good 
lattice  mismatch  with  the  Ag(lll)  surface  provided  that  a  2V2  x  2V2  r45°  cell  is  considered.  We 
also  analysed  the  interface  of  the  rutile  (100)  face  with  Ag(  111)  where,  however,  a  substantial  strain 
is  still  present.  In  both  cases  the  bonding  at  the  interface  is  very  weak  for  a  stoichiometric  oxide, 
and  the  charge  transfer  is  negligible.  Stronger  interfaces  are  formed  when  reduced  phases  of  the 
oxide  are  considered.  In  particular,  in  agreement  with  previous  work,  we  find  that  the  formation 
energy  of  a  single  O  vacancy  or  of  an  array  of  O  vacancies  has  a  lower  cost  at  the  interface  with  Ag 
than  on  the  bare  titiania  surface. 

(2)  Once  the  nature  of  the  metal/oxide  interface  has  been  determined,  we  have  considered 
the  process  of  incorporation  of  O  atoms  in  the  metal  electrode  or  of  Ag  atom  in  the  oxide  layer. 
Both  processes  imply  relatively  high  formation  energies.  This  is  due  on  one  side  to  the  relatively 
modest  affinity  of  Ag  for  O  (AH°(Ag20)  =  -31.1  KJ/mol)  compared  to  that  of  Ti  (AH°(Ti02)  = 
-944.0  KJ/mol)  and  on  the  other  side  to  the  small  size  of  the  interstitial  cavities  present  in  the  bulk 
of  anatase  and  rutile  Ti02.  The  Ag  atom  inside  these  cavities  is  compressed  and  its  steric  repulsion 
with  the  surrounding  is  significant.  The  diffusion  of  Ag  and  O  atoms/ions  across  the  interface  is 
therefore  a  thermodynamically  unfavorable  process  which  can  occur  only  at  high  temperatures  or 
under  the  effect  of  an  external  electric  field  and  that  could  be  greatly  facilitated  by  the  presence  of 
point  defects  (vacancies),  extended  defects  (grain  boundaries),  or  by  the  formation  of  an  amorphous 
phase  of  the  oxide  where  larger  cavities  may  exist. 

(3)  Assuming  that  Ag  atoms  have  penetrated  into  the  bulk  oxide,  we  have  considered  their 
electronic  nature.  We  have  studied  both  isolated  Ag  atoms  in  interstitial  or  substitutional  to  Ti 
positions,  or  higher  concentrations  of  Ag  atoms  resulting  in  the  formation  of  dense  arrays 


17 


(nanofilaments).  In  all  cases  we  observe  a  tendency  of  Ag  to  be  ionized  and  transfer  the  valence 
electron  to  the  host  crystal.  This  is  in  part  a  consequence  of  the  steric  repulsion  described  above 
since  the  transfer  of  the  Ag  5s  electron  to  the  empty  Ti  3d  states  strongly  contributes  to  reduce  the 
repulsion.  We  also  observe  a  very  efficient  screening  of  the  positively  charged  Ag  ions  so  that  they 
remain  in  this  electronic  situation  even  when  the  Ag-Ag  distances  are  short  and  extended  chains  are 
formed.  Even  when  extra  electrons  are  added  to  the  system  they  do  not  enter  in  the  Ag-Ag 
conducting  state  but  rather  prefer  to  occupy  localized  Ti  3d  states  at  the  bottom  of  the  Ti02 
conduction  band.  This  suggests  that  for  an  Ag  filament  inside  Ti02  to  exhibit  conductive  behavior  a 
higher  density  of  Ag  atoms  is  required  in  order  to  form  a  thicker  metallic  wire.  This  is  hardly 
possible  in  the  regular  bulk  lattice  of  Ti02  without  creating  a  structural  breakdown. 

Aspects  which  have  not  been  considered  in  this  work  and  that  could  significantly  affect  the 
phenomena  described  are  the  role  of  O  vacancies  and  their  accumulation  at  the  interface,  the  change 
of  the  properties  when  an  amorphous  oxide  is  considered  and  the  role  of  an  external  electric  field  in 
promoting  diffusion  of  ionic  species  across  the  interface.  All  these  effects  are  expected  to  lower 
substantially  the  energy  barriers  involved  in  ion  diffusion  and  in  the  generation  of  metal  dendrimers 
inside  the  oxide  structure. 
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